Overspill avalanching in a dense reservoir network 
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Sustainability of communities, agriculture, and industry is strongly dependent on an effective storage and 
supply of water resources. In some regions the economic growth has led to a level of water demand which 
can only be accomplished through efficient reservoir networks. Such infrastructures are not always planned at 
larger scale but rather made by farmers according to their local needs of irrigation during droughts. Based on 
extensive data from the upper Jaguaribe basin, one of the world's largest system of reservoirs, located in the 
Brazilian semiarid northeast, we reveal that surprisingly it self-organizes into a scale-free network exhibiting 
also a power-law in the distribution of the lakes and avalanches of discharges. With a new self-organized- 
criticality-type model we manage to explain the novel critical exponents. Implementing a flow model we are able 
to reproduce the measured overspill evolution providing a tool for catastrophe mitigation and future planning. 



Water is simultaneously the most relevant commodity for 
life and the primary medium through which people and the 
environment are affected by natural disasters 0J[6). Accord- 
ing to the Food and Agriculture Organization of the United 
Nations, the water demand in the last century has grown more 
than twice the population increase and as forecast, by 2050, 
about 2 billion people will be living in water scarcity [7 1. This 
is not only related to the availability of freshwater but also 
with the lack of proper storage and distribution infrastructures. 
It is thus of paramount interest to understand how to guaran- 
tee an even spatial and temporal distribution of water, avoid- 
ing the imbalance between supply and demand. The first step 
in this direction, which we address here, is to analyze how 
the interplay between the rainfall and evaporation affects the 
filling level in the reservoirs and the network. 

In semiarid environments like the Brazilian northeast it is 
common to have droughts drastically affecting the crops. The 
weather in the region is characterized by intermittent rains and 
long periods of water scarcity. The water supply, in the dry 
season, is assured by a dense set of surface reservoirs, which 
are essential for the sustainability of the region, since during 
the rainy season water is stored to be distributed afterwards. 
A proper reservoir size is crucial: while smaller reservoirs are 
typically unable to supply enough water throughout the entire 
dry season, due to the low storing capacity and high losses 
through evaporation and infiltration, too large reservoirs can 
have serious consequences for life sustainability. For exam- 
ple, if the storage capacity in a basin exceeds the volume 
equivalent to three times the average annual runoff, its water 
will not be sufficiently renewed |8|. Not only the properties 
of each reservoir are relevant but also the way they are inter- 
connected. For example, upstream reservoirs may retain a sig- 
nificant part of the collected water which enables an efficient 
spatial distribution by ensuring a continuous connection with 
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other reservoirs even during dry seasons, when usually no wa- 
ter flows through the river network [9-11]. The available data 
from interferometric radar and weather stations lfT2llT5l pro- 
vide the necessary ingredients for a systematic analysis of the 
existing network, which is fundamental to further planning for 
improvement. We introduce here a hydrological model for the 
evolution of water resources in the network and apply it to 
the upper Jaguaribe basin, a semiarid region in the Brazilian 
northeast, for which we have collected information about the 
rainfall and the distribution of reservoirs over two decades. 



I. THE RESERVOIR NETWORK 

Located in the state of Ceara, northeast of Brazil, the upper 
Jaguaribe basin has a catchment area of about 25 000 km 2 
and almost 4 000 reservoirs (see Fig. [TJ. Characterized by 
a broad reservoir-size distribution, ranging from 2 500 m 3 
to 2 billion m 3 , this basin stores water for more than half- 
million people. It has an average annual rainfall of 860 mm 
and a potential evaporation as high as 2 000 mm. The 
weather conditions together with shallow soils above crys- 
talline bedrocks make the rivers intermittent and leads to 
runoff rates below 7% of the precipitation. The main river is 
the Jaguaribe, which contributes to the Oros reservoir located 
at the basin outlet, the black reservoir in Fig. [T] (northeast). 
The latter, with a storage capacity of 1 940 million m 3 , is the 
second largest reservoir of the state and an important water 
reserve, which yields about 630 million m 3 /year with an an- 
nual reliability level of 90%, i.e., 90% probability of being 
able to supply this amount every year. 

Analyzing satellite images from the wet years (see section 
Materials and Methods), we have identified the geographical 
location and size of each reservoir, as well as the connection 
between them. The location of the full set of reservoirs in the 
Jaguaribe basin is represented in Fig.[T]C. The color scheme 
stands for the distance from the outlet (Oros), in units of reser- 
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voirs. As shown in Fig.[T]D, to characterize the network we 
schematize it as a graph where reservoirs are nodes and the 
connections between them are links. The links are directional 
since flow goes from upstream to downstream, based on the 
height difference between reservoirs. The collected data dis- 
close a network of reservoirs organized in a tree with low 
branching rate. All reservoirs have a single outlet connection 
but with a broad distribution of inlets. While most upstream 
reservoirs have no inlet connection, for the ones downstream 
typically multiple connections exist (see Fig.[T]D). The degree 
k of the reservoir (i.e., number of connections) varies from one 
to around 400 and the degree distribution Pk(k), is a power 
law, Pt(k) ~ k~V, with a degree exponent y = 2.3 ± 0.2, as 
shown in the inset of Fig. [2] 

The satellite images from the wet years allow to measure 
the area (A) of each reservoir disclosing a power-law distribu- 
tion of areas, P A (A) ~A~ X , with X = 2.00 ±0.05 (see Fig.^. 
From the area, the storage capacity (V) can also be estimated. 
Molle [ 3 1 1 has studied 417 reservoirs of the Brazilian semiarid 
region, a 970 000 km 2 domain in the northeast and southeast 
of Brazil, and proposed an area/volume relation given by, 

A = Cd {d) ' (1) 

where c and d have been established empirically to be c = 2.7 
and d = 1500. For the upper Jaguaribe basin, Mamede (9) 
tested this equation for 21 reservoirs showing satisfactory 
agreement. Since the capacity of a reservoir cannot be di- 
rectly measured with the satellite imagery, we apply the same 
relation (Eq.[T]i to the remaining reservoirs. 

H. HYDROLOGICAL MODEL 

The planning of an efficient reservoir network requires the 
analysis of its response to typical weather conditions. The 
amount of water in a reservoir is a dynamical quantity, with a 
positive contribution from precipitation and upstream inflow, 
while mechanisms such as groundwater infiltration, evapora- 
tion, and overflow tend to reduce it. To quantify the precipita- 
tion and evaporation, meteorological time series for monthly 
evaporation level and rainfall (volume per area) have been ob- 
tained for the 131 weather stations of the region [16|. From 
the topology of the basin, the watershed lines have been traced 
and the subbasin of each reservoir determined. The daily av- 
erage of rainfall and evaporation in the subbasin is determined 
from the closest station. The net evaporation is estimated 
based on the mean reservoir area, while the contribution of 
the precipitation to the reservoir level accounts for both direct 
rainfall in its area and a fraction of the one falling in the area 
of the subbasin drained into it. The latter fraction is denoted as 
runoff and depends on factors such as soil constituents, vege- 
tation coverage, rain intensity (both actual and preceding), and 
topography fTTl - EOl . For simplicity, we take the same value of 
the runoff (0.04) for the entire region, i.e., 4% of the rainfall 
in the subbasin, outside the reservoir, is considered to drain 
into the reservoir. As discussed below, a good agreement is 



found for the time evolution of the water volume in different 
reservoirs with this uniform runoff coefficient. 

To model the time evolution of the water resources in the 
reservoir network we consider that initially each reservoir is 
at 20% of its storage capacity (presented results are not af- 
fected by the initial conditions). Following the time series for 
the rainfall and evaporation, the inflow of water to each reser- 
voir is computed on a daily basis accounting for three con- 
tributions: direct rainfall, water drained from the subbasin, 
and inflow from reservoirs upstream. Reservoirs are updated 
from upstream to downstream and the overflow is obtained 
from the difference between the total amount of water and the 
storage capacity. The exceeding water is transferred to the 
reservoirs downstream, where we neglect transmission losses 
through the river transport. When an overflow occurs, water 
is transferred to the next reservoir which can trigger a cas- 
cade of overflows of other reservoirs downstream. We denote 
as avalanche this type of events and we measure their size as 
the number of connected reservoirs which overflow within the 
same day. 



III. WATER TRANSPORT 

For simplicity, let us start by considering a uniform distri- 
bution of rainfall and evaporation over the entire region during 
two decades. At the subbasin level, a competition takes place 
between the increase of water due to the rain and upstream 
inflow and the evaporation. Reservoirs where the influx is 
larger than the outfiux, after some initial filling time, are full 
and constantly overflowing, i.e., transferring water to the one 
downstream. On the other hand, reservoirs where losses are 
not compensated by incoming water, essentially dry out and 
overflow never occurs. As an example, we take a level of 
precipitation of 1.825 mm and an evaporation of 6.832 mm 
(per unit area and day), corresponding to the daily averages 
for the period between 1991 and 2010. When we initialize 
the system with all reservoirs at 20% of their capacity, after 
15 years, every day, one avalanche affecting 1258 reservoirs 
occurs. This scenario is not observed in reality since both the 
precipitation and evaporation are asymmetrically distributed 
in space and time. Therefore, a meaningful modeling of the 
water transport requires a proper quantitative description of 
these two key processes. 

Inserting the time series for rainfall and evaporation from 
1991 till the end of 2010, we have followed the dynamics of 
the reservoir network and measured the size of the avalanches 
as well as the volume of water storage at each reservoir. We 
define the avalanche size s as the number of connected reser- 
voirs from which water overflows during the same day. Over 
these two decades, in the rainy seasons, many avalanches are 
observed with sizes ranging over three orders of magnitude. 
Most of the avalanches solely affect two reservoirs, but large 
ones also occur (see Supporting Video). The avalanche-size 
histogram H(s), shown in Fig. [3] scales with the avalanche 
size as H(s) ~ s~ a , with a = 1.9 ±0.2. The power law be- 
havior is related to the absence of a characteristic scale in the 
avalanches (see section Stochastic Model). In the considered 
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time interval, the largest observed avalanche occurred during 
the day 47th of 2004 and involved about 2000 reservoirs, i.e., 
more than 50% of the entire set. In the inset of Fig. [3] we 
also show the distribution of the area P^, affected during an 
avalanche. The affected area Af is defined as the total area of 
catchments related to the overflowing reservoirs. This distri- 
bution scales as P Af ~ AJ P , with /3 = 0.56 ± 0.08. This low 
value of j3 reveals that large scale events evolving extended 
catchment areas also occur. 

A constant runoff of 0.04 has been considered. The expo- 
nent of the avalanche-size histogram is resilient over a wide 
range of runoff coefficients (below 0.05). For large values of 
runoff some reservoirs have a large inflow of water when com- 
pared to the outflow and, once full, they overflow every day. 
While the exponent of the avalanche-size histogram does not 
depend on the runoff coefficient, to reproduce the measured 
water volume of each reservoir, a proper value is required. As 
pictured in Fig. [4] with the considered runoff, the water vol- 
ume is satisfactorily reproduced over most of the considered 
time range. Slight differences between the simulation and real 
values occur after long periods of intensive rain. As discussed 
by Koster et al. 1 19|, in this case the runoff coefficient could 
have been affected by the antecedent rainfall. As observed for 
2004, after long periods of intensive rain, the wet soils en- 
hance the runoff, promoting the drainage of a larger volume 
of water from the subbasin to the reservoir and thus larger 
avalanches. 



IV. STOCHASTIC MODEL 

The most striking result unveiled by our hydrological model 
is the power-law avalanche-size distribution. This result re- 
sembles the fingerprint of self-organized critical (SOC) sys- 
tems disclosed by P. Bak [22 1 though with a different exponent 
than SOC on scale-free trees l23l |24l . SOC is characterized 
by branching, where the discharge into neighbors might trig- 
ger also their discharge under slow driving. However, none 
of these features appears in the reservoir networks. Rain falls 
synchronously on several reservoirs and, while the inflow con- 
nections are scale free, the outflow occurs through one single 
connection (no branching). Furthermore, water coming from 
the synchronous overspill at two (or more) upstream reser- 
voirs joins on downstream ones and increases the amount of 
water transported through a single branch. 

Let us consider the following model for the same set of 
reservoirs from the upper Jaguaribe basin. Each reservoir i 
is characterized by its amount of water m, . Iterative ly, a unit 
volume of water is added to randomly selected reservoirs. 
Though on average each reservoir receives the same amount 
of water, at iteration i a reservoir can be selected more than 
once. To each reservoir a threshold f, is assigned such that if 
the amount of water in the reservoir exceeds this threshold, 
fiij > tj, overspill occurs and a fraction 1 — / of water is trans- 
ferred to the downstream reservoir while the remaining water 
/ is dissipated (evaporation and transmission losses). Since 
the inset of Fig.[5]shows a linear relation between the capacity 
of a reservoir and its number of connections (degree), we take 



ti as equal to the degree. In the main plot of Fig. [5] we show 
the avalanche-size distributions for different values of dissi- 
pation / obtained simulating the above model on the 3 978 
reservoirs identified in the upper Jaguaribe basin. For a wide 
range off, results are consistent with a power-law distribution 
extended over three orders of magnitude with the same expo- 
nent a = 1.9 obtained with the hydrological model in Fig. [3] 
(solid line). 



V. FINAL REMARKS 

In summary, we found some unexpected scale-free behav- 
ior in the dynamics of water transport on a reservoir network, 
using available data for the geographical distribution of reser- 
voirs as well as the time series of rainfall and water evapora- 
tion. The hydrological model for the in- and outflow balance 
at each reservoir reveals an excellent qualitative and, in most 
of the cases, quantitative agreement for the upper Jaguaribe 
basin over two decades. The analysis of these reservoirs, with 
satellite imagery, discloses a tree-like scale-free network with 
a power-law distribution of capacities. Our model for the wa- 
ter transport shows a cascade dynamics of connected spilling 
over reservoirs characterized by overflows extended over three 
orders of magnitude in length scales. The introduced stochas- 
tic model grasps the main features observed with the hydro- 
logical model. With this scheme it is now possible to charac- 
terize the role of the network topology, spatial and time cor- 
relation of the rainfall, and evaporation, on water flow. The 
proposed tools aim also to help the planning of new reservoirs 
by numerically testing their impact under typical weather con- 
ditions, avoiding either scarcity over dry seasons or floods in 
rainy ones. Despite the degree of agreement between numer- 
ical and real data, further improvement of the model may re- 
produce even more details. Future work can account for the 
temporal evolution of the network itself and its impact in the 
region [25]. In this study the water transport between reser- 
voirs has been considered to occur within one day and without 
transmission losses, a more detailed description of the mass 
transport could also be included [26 , 27 1. Upstream reservoirs 
also serve as sediment detention basins, retaining a consider- 
able amount of sediment generated within the catchment and 
extending the lifetime of larger reservoirs located downstream 
EQT), affecting the storage capacity and, consequently, the 
entire dynamics l28l l29l . Understanding the impact of such 
mechanism may also be an interesting focus of research. 



VI. MATERIALS AND METHODS 

To identify the reservoirs in the network and quantify their 
maximum capacity we analyzed images from three wet years, 
namely, 2004, 2008, and 2009. These images, collected by 
the satellites LANDSAT (1, 2, 3, 5, and 7) and CBERS (2 and 
2B), were obtained from the website of the "Instituto Nacional 
de Pesquisas Espaciais" ll30l . The LANDSAT generates im- 
ages in seven bandwidths: three in the visible region, one in 
the near infrared, two in the mid-infrared, one in the thermal 
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infrared, and a panchromatic one. By properly selecting refer- 
ence points, the obtained images have been georeferenced and 
prepared to be used in the Geographical Information System 
(GIS). Through unsupervised classification techniques, which 
do not require any further information besides the images, the 
data have been processed and pixels grouped by spectral fea- 
tures to identify the ones corresponding to water reservoirs. 
An additional handmade filtering was necessary to avoid mis- 
leading information due to the presence of clouds. By super- 
posing the images for different years, it was possible to lo- 
cate each reservoir and trace its perimeter, obtaining the max- 
imum area. With the relation proposed by Molle [31 1, Eq. [T] 
the maximum capacity of each reservoir has been estimated. 
The drainage network was obtained from the analysis of high- 
resolution elevation data, collected by the Shuttle Radar To- 
pography Mission (SRTM) from NASA HUE], by tracing 
the watershed lines, and afterward each reservoir was located 
in this network. The dataset with the network and geograph- 



ical localization of each reservoir is available as Supporting 
Information. The daily average of rainfall for each reservoir 
is obtained from the closest weather stations and the evap- 
oration is estimated from the monthly average over several 
decades fl6l . The pan evaporation data accounts for the evap- 
oration on the surface of the reservoir (about 70% of class A 
pan evaporation) ||34l , for the infiltration, as well as for the 
water extracted for several uses (about 30%) l8ll3~Tl. 
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FIG. 1. Geographical location of the Oros Basin, in the Brazilian northeast (A). Located in the state of Ceara (B), the Upper Jaguaribe 
Basin consists of 3 978 water reservoirs with a total catchment area of about 25 000 km 2 . The reservoirs have been identified through satellite 
imagery from three wet years, namely, 2004, 2008, and 2009. In (C) we mark the reservoirs where the color scheme represents the distance to 
the outlet in units of intermediate reservoirs. The network can be schematized by a graph (D) where nodes represent the reservoirs and links 
stand for the connections between them. For the sake of graphical clearness, reservoirs without incoming links have not been included. 
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FIG. 2. The reservoirs in the upper Jaguaribe basin are characterized 
by a broad distribution of areas and connections. Main plot: The 
distribution of areas Pa (A) scales according to a power law, Pa (A) ~ 
A~ x , with X = 2.00 ±0.05. Inset: The degree distribution P k {k) of 
the current reservoir network follows a power law, Pk(k) ~ k~^, with 
7 = 2.3 ± 0.2. The size and location of the reservoirs were measured 
from satellite images (see section Materials and Methods). 
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FIG. 3. Power laws over more than three orders of magnitude, 
revealing a lack of a characteristic size on the avalanches and af- 
fected areas. Main plot: The avalanche-size histogram scales as 
H(s) ~ s~ a , with a = 1.9±0.2, where the avalanche size s is defined 
as the number of connected reservoirs with overflow. Inset: The dis- 
tribution of affected area P Af (A /) ~ AJ P , with j8 = 0.56 ± 0.08. The 
area A t corresponds to the total area of catchments related to the 
avalanche. Both results have been numerically obtained for a runoff 
coefficient of 0.04 and with the real time series for rainfall and evap- 
oration, between 1991 and the end of 2010. 
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FIG. 4. Model comparison with empirical results for the temporal 
evolution of the water volume at the Oros reservoir, revealing a good 
quantitative agreement between both. Numerical results have been 
obtained with simulations for a runoff coefficient of 0.04 and the real 
data corresponds to the time series of the water volume in the interval 
between 1991 and 2011 1211 . The larger differences are observed 
after long periods of intensive rain (like the beginning of 2004) since 
the wet soils improve the amount of water drained from the subbasin 
to the reservoir. The differences can also be due to the construction 
of new reservoirs during this two decades. 




FIG. 5. Power laws with the same exponent than the ones for the hy- 
drological model are observed with a simplified model on the reser- 
voir network from the upper laguaribe basin. Main plot: Avalanche- 
size distribution obtained for different values of dissipation, namely, 
0.4, 0.5, 0.6, and 0.7. Results have been averaged over 10 2 samples 
of 10 s iterations each. The solid line stands for a power law with 
exponent 1.9 as in the main plot of Fig. [2] Inset: Dependence of 
the average reservoir capacity on the degree of the reservoir, where 
the solid line corresponds to a linear relation. The capacity of each 
reservoir have been estimated from satellite imagery (see section The 
Reservoir Network). 



